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We employ effective mass theory for degenerate hole-bands 
to calculate the acceptor binding energies for Be, Mg, Zn, Ca, 
C and Si substitutional acceptors in GaN and A1N. The calcu- 
lations are performed through the 6x6 Rashba-Sheka-Pikus 
and the Luttinger-Kohn matrix Hamiltonians for wurtzite 
(WZ) and zincblende (ZB) crystal phases, respectively. An 
analytic representation for the acceptor pseudopotential is 
used to introduce the specific nature of the impurity atoms. 
The energy shift due to polaron effects is also considered in 
this approach. The ionization energy estimates are in very 
good agreement with those reported experimentally in WZ- 
GaN. The binding energies for ZB-GaN acceptors are all pre- 
dicted to be shallower than the corresponding impurities in 
the WZ phase. The binding energy dependence upon the 
crystal field splitting in WZ-GaN is analyzed. Ionization lev- 
els in A1N are found to have similar 'shallow' values to those 
in GaN, but with some important differences, which depend 
on the band structure parameterizations, especially the value 
of crystal field splitting used. 



I. INTRODUCTION 

Wide-band gap III-V nitrides, particularly Ga-, Al- 
and InN, and their semiconductor alloys, are materi- 
als currently under intense study. Some of their most 
promising applications in optoelectronics devices are 
for instance the fabrication of blue/green LEjD's,El laser 
diodesp and 'solar-blind' UV photodetectors.B The per- 
formance improvements of these and related optoelec- 
tronic devices depend strongly on the features of the in- 
trinsic and extrinsic impurity defects in the nitride com- 
pounds. For example, defects and impurities provide free 
carriers under suitable conditions. Therefore, knowing 
the accurate position of the donor and acceptor levels 
of these systems is an issue of great importance in the 
understanding of optical properties and practical appli- 
cations of these nitrides. 

At present, Mg and Zn are the impurity materials most 
widely employed in the p-doping of GaN. The experimen- 
tal thermal ionization energy (acceptor binding energy) 
associated with Mg is estimated at 250 meV.Q The high- 
est doping achieved reaches hole concentrations af ap- 
proximately 3 x 10 18 cm" 3 at room temperature.El It is 
also known that in order to activate the dopants and 
improve the p-type conductivity, the samples must be 



treated with low energy electron beam irradiation, fut. 
nace annealing, or rapid thermal annealing after growth.tl 
On the other hand, Zn doping seems to be inefficient b 
cause of its relatively deep ionization energy (340 meV) 
Other dopants have been considered, but experimen- 
tal problems like instability and/or hole compensation 
due to the formation of acceptor-H neutral complexes 
is still at issue. Estimates for the binding energies of 
several substitutional acceptors in GaN have been ob- 
tained in the_past, mostly through photoluminescence 
(PL) spectra.u However residual impurities and defects 
in this material complicates the identification of these 
levels. In contrast, little is known about the doping and 
spectrum of impurity levels in A1N. In fact, no conclu- 
sive results for the doping of A1N with sufficiently high 
conductivity have yet been reported. 

Apart from the question of successful p-doping in GaN 
and A1N using various impurities, there are still at least 
two other important issues that are under scrutiny. The 
first one is related to the determination of the origin of 
the chemical shift observed in the acceptor spectrum lev- 
els in GaN, apparently induced by the differences in the 
cores of the various impurity atoms, and some possible 
lattice relaxation around the impurity atom.Q The sec- 
ond question is whether acceptors with smaller binding 
energies (< 230 meV) exist for wurtzite and zincblende 
GaN. The occurrence of relatively large ionization ener- 
gies for acceptors in GaN has been attributed in part 
to the fact that the III-V nitrides are more ionic than 
other III-V compounds (such as GaAs, GaP and InP), 
for which the acceptor binding energies are an order of 
magnitude smaller than those found in GaN. It has also 
been suggested that the enhanced binding energies found 
for some acceptors like Zn and Cd_ are associated to the 
relaxation of the d-electron core.Ll On the other hand, 
impurities without d-electron states, such as Mg, C, and 
Si, appear to induce rather shallow .acceptor levels. In- 
deed, very recently, Park and Chadicl examined the sta- 
bility of acceptor centers in GaN, A1N and BN using first 
principles calculations. They concluded that the small 
bondlengths in III-V nitrides inhibit large lattice strain 
relaxations around impurities (mainly Be, Mg and C), 
giving rise to relatively shallow states for these species. 
This would suggest that a similar lack of relaxation ac- 
companies other substitutional impurities in these hosts, 
producing relatively shallow levels, as long as there are 
no d-cores close to the valence band energies. 

Very recently, the formation energies and impurity lev- 
els for a few donor and acceptor species have been studied 
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theoretically by several groupsJUll employing quantum 
molecular dynamics schemes and total energy calcula- 
tions in the local density approximation of density func- 
tional theory. In general, consistency is found among 
those groups, as well as with experimental reports for 
some impurity levels, such as Mgo a acceptors (Xy in- 
dicates the ion X substituting in the Y site). However, 
this is not the case for other acceptors like Cjv, where 
discrepancies of factors of three exist among theoretical 
values. Although the calculated energy levels for these 
approaches appear reliable for most cases, the impurity 
levels reported for some acceptors are close to the sys- 
tematic error bars introduced in the calculations. The 
delicate and complex nature of these calculations, which 
require intensive computations, suggest that alternative 
methods should be explored in the study of impurity lev- 
els in these systems. There is also, no doubt, the need 
for new careful experiments in the better-characterized 
materials now available, to clarify these discrepancies. 
The features of the acceptor states in the different crys- 
tal phases, wurtzite (WZ) and zincblende (ZB), have not 
been discussed either. In order to address these ques- 
tions, we present here a contribution towards the the- 
oretical treatment of the impurity levels in GaN and 
A1N based in the effective mass approach for degenerate 
bands. 

In this paper we report effective mass theory calcula- 
tions of the acceptor binding energies for various impu- 
rity atoms in GaN and A1N for both crystal structures, 
WZ and ZB. Particular attention has been paid to chem- 
ical shifts introduced by the foreign atoms. An acceptor- 
pseudopotential model is used to take into account this 
effect. The approach used here is based on the effective 
mass theory (EMT) for degenerate bands. Well param- 
eterized valence band structure calculations are used as 
input. The results obtained, with no adjustable param- 
eters, are in very good agreement with experiments, as 
we will see below. Inevitably, the application of even a 
simple hydrogen-like model of acceptor states in group 
III-V semiconductors is more complicated than for ide- 
alized semiconductors with a single, isotropic and spin- 
degenerate valence band. The complications are due in 
part to the band warping and sixfold degeneracy or near- 
degeneracy of the valence band structure close to the T 
point (k = 0). Since the perturbing potential introduced 
by the foreign atoms can be seen to zeroth order as pure 
Coulomb-like, the problem can be seen as a generalized 
hydrogenic problem, where the kinetic energy of a hole, 
in the rather complicated valence band structure of the 
III-V materials, is properly described by a 6 x 6 matrix 
Hamiltonian which describes well the dispersion features 
of the various hole bands. The EMT calculations of the 
binding energies of Be, Mg, Zn, Ca, and C acceptor im- 
purities are shown to be in very good agreement with 
the available experiments, and consistent in general with 
other theoretical calculations employing other methods 
(with the exceptions discussed above for C, for example). 
The applicability of EMT for the calculation of impurity 



levels with 0.2-0.4 eV binding energies is then verified 
post facto, likely due to the large bandgap in these mate- 
rials, which yields negligible mixing of conduction band 
states. 

Additionally, we find that the binding energies for ac- 
ceptors in the ZB structures are predicted to be shallower 
than their counterparts in the WZ structures, suggest- 
ing that doping of ZB material would be of significant 
practical advantage. We notice that the difference on 
parameters, mainly the existence of a crystal field split- 
ting for the WZ nitrides, strongly affects the band mix- 
ing and correspondingly the binding energies in the two 
polytypes. It is also likely that differences in the hole 
masses contribute to the calculated different binding en- 
ergies. Although substitutional impurities do not repre- 
sent a strict test of the different band parameterizations, 
the subtle interplay of the different valence bands on the 
resulting binding energies provides an interesting overall 
consistency check of the parameterized band structure. 

The paper is organized as follows. In section II, we 
present the characteristics of the generalized acceptor 
problem. The explicit matrix form of the ZB and WZ va- 
lence band Hamiltonians are also given there. The trial 
form of the envelope wave functions is presented in sec- 
tion III. The impurity pseudopotential model is discussed 
in detail in section IV. The correction due to polaron ef- 
fects is briefly described in section V. The results and 
discussion are given in section VI, and the conclusions in 
section VII. 



II. GENERALIZED SHALLOW ACCEPTOR 
PROBLEM 

Substitutional impurities with one fewer valence elec- 
tron than the host atom of the pure crystal introduce 
well localized acceptor states lying just above the top of 
the valence band structure. The theory of shallow donor 
and acceptor states in semiconductors has been reviewed 
in detail by Pantelideso We assume, as usual, that all 
acceptor levels in the semiconductor are described within 
the effective mass theory for degenerate band structures 
by the following matrix equation 



#(r)F(r) 



[H{v) + t/(r)l]F(r) 
EF(r) , 



(1) 



where H (r) is the full acceptor Hamiltonian with eigen- 
values E for the acceptor states. Here, 7i(r) is the Hamil- 
tonian properly constructed from crystal symmetry con- 
siderations which entirely describes the spectrum and 
eigenvalues of a hole near the valence band extrcmiun. 
at the r point. Symmetry invariance group theor\ i. lii 
and kp perturbation theory for degenerate bandaL3'E2 
has been used to derive the proper effective-mass Hamil- 
tonian for strained semiconductors depending upon the 
crystal structure. 
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The potential U(r) is the perturbation produced by 
the acceptor-ion on the otherwise pure and periodic host 
crystal. In a simple idealized case, U(r) is taken to be the 
Coulomb potential U(r) — e 2 /e |r|, where e is the static 
dielectric constant of the crystal, e(q = 0, u> = 0), rep- 
resenting a point charge in a dielectric medium. Notice 
that the screening of the simple hydrogenic potential by 
a dielectric function e(q) has been considered in the past 
as an approach to consider the contribution to the accep- 
tor spectpum of the short range potential from the real 
impurity.ta Although this model gives an insight into the 
specific character of the different atomic acceptor levels, 
the model results in a generic value for all the impurity 
defects. This, clearly, neglects the chemical signature 
of the foreign atoms in the. host material (the so-called 
central-cell contribution). o Given these limitations, we 
employ instead an ab initio pseudopotential U ps (r) corre- 
sponding to the difference between the bare model poten- 
tial of the impurity and the host atoms. Since the chem- 
ical correction induced by different species is expected to 
be small and because the pseudopotential used is fairly 
smooth and without discontinuities, the effective mass 
approach is expected to yield an appropriate description 
of the system. More details on the impurity potentials 
used are given in section IV. 

In Eq. ([jj), 1 is the 6x6 unit matrix and F(r) is a col- 
umn vector whose Fj (r) elements characterize the enve- 
lope function which modulates the Bloch functions <pj (r) 
of the unperturbed crystal at the top (k ss 0) of the va- 
lence structure. Correspondingly, the wave functions for 
the shallow states are given by 

6 

tf(r)=X>i(r)0i(r)- ( 2 ) 
i=i 

The trial form chosen for the envelope functions Fj(r) 
is discussed in detail in section III. In the following sub- 
sections we describe briefly the explicit form of the hole 
Hamiltonian H(r) for the two crystal polytypes (WZ and 
ZB), in which the bulk GaN and A1N semiconductors 
grow. 

A. Wurtzite valence band Hamiltonian 

In order to consider the motion of a carrier at the top 
of the valence band in a wurtzite semiconductor we must 
take into account its six-fold rotational symmetry, which 
induces a crystal field splitting. Moreover, in the case 
of spin-orbit interaction, the Ti5 level splits into the Tg 
state, upper V*? level, and lower I\ level, corresponding 
to the heavy hole, light hole and split-off hole bands.t£l 
The appropriate effective mass Hamiltonian that reflects 
those features of the WZ GaN bulk crystal should be 
described-thus by the Rashba-Sheka-Pikus Hamiltauian 
(RSP),E_§Ej as discussed recently by Sirenko et aZ.EJ In 
the vicinity of the valence band maximum, and to second 



order in k, the six states (including the spin index) of the 
RSP Hamiltonian for unstrained WZ structures can be 
explicitly written in a matrix representation as follows: 
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where 



F = 


X + + Ai + A 2 


G = 


\ + + A 1 -A 2 


A = 


A 1 k 2 z + A 2 k 2 ± 


= 




H = 


i{Afik z k + + A?k+ 


I = 


i{A§k z k+ — A-jk + 


K = 


A 5 k 2 + 


A = 





with k\ — k\ + k 2 , and k± — k x ± ik y . Here, Ai 
corresponds to the energy splitting produced by the 
anisotropy of the hexagonal symmetry, A 2 = Ago /3 and 

A3 = Aio /3 are the energy splittings for the z and per- 
pendicular jdirections produced by the spin-orbit (SO) 
interaction.^ The A constants are related to the inverse 
of the hole masses, in units of fi 2 /2m , where m is the 
bare electron mass. Notice that when the linear terms in 
(^) are negligible (A 7 — 0; which is in fact nearly the case 
in GaN and A1N), the RSP Hamiltonian has complete 
inversion symmetry. This symmetry allows for helpful 
simplifications in dealing with the acceptor problem in 
the envelope function framework, as we discuss below. 

B. Zincblende valence band Hamiltonian 

In the case of semiconductors with the ZB structure, 
the hole wave functions characterizing the sixfold degen- 
erate Ti5 state split, due to the effects of spin-orbit inter- 
action, into the fourfold degenerate Tg states correspond- 
ing to the heavy apd light hole bands, and the spin-split 
off hole states TjS3 The Hamiltonian which takes into 
account all these features of the cubic symmetry for ZB 
semiconductacs is the well known Luttinger-Kohn Hamil- 
tonian (LK),E3 which at the valence band extremum, and 
to second order in k, is expressed in terms of only four 
empirical parameters — the so-called Luttinger-Kohn pa- 
rameters 71,72 and 73, and the spin-orbit splitting A Q . 
Thus, the LK Hamiltonian TLzb is written in matrix form 
as follows 
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(5) 

where 

L = —2\fii^ 2 k z k- 
M = V3 l3 (k 2 x - k 2 y ) - 2V3ij 3 k x k y 

P = lik 2 - 72 (2k 2 z - k\) 
Q = ll k 2 + l2 (2k 2 z -kl) 

R=-^i{P-Q) 
S = -iV2M 
W=^(2P + Q) + A 

Y = ±(P + 2Q). (6) 

Here, L, M, P and Q are in units of h 2 /2m . Notice 
that the higher symmetry of the ZB structure produces 
a much simpler H(r) and fewer parameters than for the 
RSP case. In both cases, the operator Tt(r) is obtained 
via the usual transformation k a — > igf- in Ti(k). 

III. TRIAL FORM FOR THE ENVELOPE 
FUNCTIONS 

To solve the effective mass equation for degenerate 
bands, Eq. ([!]), we use the fact that the effective mass 
Hamiltonian is invariant under inversion with respect to 
the origin, so that the envelope functions Fj(r) can be 
chosen to have definite parity. Since the features of the 
acceptor problem are rather like those of a hydrogenic- 
like problem, it has proved convenient to choose the en- 
velope functions basically as an expansion in spherical 
harmonics and a linear combination of hydrogcnic-likc 
radial functions. In particular, we have chosen the fol- 
lowing explicit form: 

summing over all I even (or odd), and with radial func- 
tions for a given hole band j and angular momentum 
quantum number I of the form 

N 

//(r) = ^i/Ve-^. (8) 

i=l 



In this work, however, we are mostly interested on the 
ground state (the highest binding acceptor state), and 
in such state only even I will contribute to the expan- 
sion — as one would expect a ground state with even 
parity. This convenient simplification can be relaxed 
straightforwardly if desired, with little effect on the re- 
sults. For numerical convenience, we find it useful to 
minimize or evaluate the acceptor binding energy choos- 
ing a-s in the progression au = aie^ -1 ), such that (3 = 
(N — 1) _1 log(ajv/ai), and the end point conditions are 
chosen as a\ — 1.2 x 10 -2 a* , and a at = 3.5 x 10 2 a* . 
Here, a* = "fie a is the effective Bohr radius, and 71 is 
defined by 

-J -(2m /h 2 )(A 2 + A 4 ) for WZ , . 

71 \ 7l for ZB ' 1 > 

such that the effective Rydberg energy is defined as E* = 
m e 4: /2h 2 "fie 2 — e 2 /2a jie . The range of ai values was 
designed to cover a wide spectrum of length scales. In the 
limit case of 71 = e = 1 (with N=25 and for I = 0,2), 
one obtains the hydrogen spectrum, so that for the first 
five states we obtain (in Rydbergs) E x = 1.0000, E 2 = 
0.2500, E 3 = 0.1111, E A = 0.0625 and E 5 = 0.0399, as 
expected. 

IV. IMPURITY ATOM PSEUDOPOTENTIAL 

As mentioned above, a simple hydrogenic (scaled 
Coulomb) potential would not yield the observed vari- 
ations in the binding energy of acceptor states for differ- 
ent impurity atom species. Photolumincscence measure- 
ments show indeed important differences in the acceptor 
binding energies in WZ GaN for different impuritiesa. To 
study those 'chemical' shifts one needs to use impurity 
potentials properly constructed to insure that their phys- 
ical properties reflect the expected shifts. The impurity 
potential here is obtained from an analytical represen- 
tation of the pseudopotential for the bare impuritj: and 
host atoms. The analytic form follows Lam et aL,c3 who 
fit the first-principles pseudopotentials developed earlier 
by Zunger et alEi in a density functional formalism. No- 
tice then that the acceptor potential is truly an impu- 
rity pseudopotential, having its origin in ab initio cal- 
culations. i-The pseudopotential for a bare atom can be 
written asEil 

[7ps(r)= ^^(r)P i -^, (10) 
1 

with 

where V^ s (r) represents the atomic core pseudopotential. 
Pi is the projection operator which picks out the compo- 
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nent of the wave function with angular momentum num- 
ber I. The constants C[, C\ and C3 are the fitted param- 
eters, with Z c and Z v representing the core and valence 
electron charges, respectively, as defined by Lam.Eil The 
first term in ( [11] ) represents a potential barrier which re- 
places the kinetic energy of the true valence states, while 
the second term arises from electrostatic screening of the 
nucleus by the core electrons and exchange-correlation 
forces. Using these pseudopotentials, the impurity model 
potential is constructed as follows. 

When the substitutional impurity atom replaces the 
host atom in the crystal, the impurity potential is de- 
fined as the difference between the impurity and host ion 
pseudopotentials. If I — 0, for instance, 

e 2 AZ e 2 

[/(r) = -A7;(r)-^, (12) 

with 

AV p ° s (r) = ±(V° s>host (r) - V p ° Stimp (r)) for Z host >Z imp . 

(13) 

Here, AZ V = z^ ost -Z^ (= 1 for single acceptors), and 
e is the dielectric constant of the host lattice. Clearly 
the first term in U(r) corresponds to the net potential 
produced by the difference between the bare core poten- 
tials of the impurity and the host; it is the short-range 
part. The last term is the long-range Coulombic poten- 
tial due to the difference in the valence charge AZ V . The 
static dielectric constant e a is introduced here to reflect 
the effect of the lattice polarizability (screening) of the 
host crystal. Notice that in this approach the net effect 
of the redistribution of charge near the impurity defect 
and the accompanying screening of the foreign charge at 
'large' distances (several lattice units) is considered fully 
in the pseudopotential definition. _ 

In a different approach, frequent in the literature,CJ the 
role of the pseudopotential is partly simulated using a q- 
dependent screening function e(q) (— » e a for q — > 0) in 
the simple hydrogenic-style impurity potential. We avoid 
using e(q) thanks to the impurity-specific pseudopoten- 
tial. We believe our approach to be better in this prob- 
lem, as it requires no further adjustable parameters and 
yields the expected chemical shifts quite accurately. 

To provide a simple and independent test of the model 
we have calculated the binding energies for several accep- 
tors in the well characterized semiconductor GaAs. The 
results are shown in Table I. The theoretical binding en- 
ergies are in excellent agreement with the experimental 
values, with no additional parameters. 

V. POLARON CORRECTION 

We should also notice that since the nitride semicon- 
ductors (GaN and A1N) are polar materials, one would 
expect that the electron-LO phonon coupling would in- 
troduce corrections to the bound states. In order to 



obtain an estimate of such correction, we assume that 
the polaron contribution to the acceptor binding energy 
close to the T point is diagonal in band index. There- 
fore the acceptor binding energies will be enhanced by 
(1 + aF(m*)/6)E*j, up to first order in the Frohlich 
coupling constant af i each hole band. This coupling 
constant is defined byE3 

wK) = (—--) fir— V' 2 > (») 

where E is one Rydberg, E* a is the ground state energy 
of the impurity acceptor without the polaron correction, 
Tilo is the LO phonon energy, and rrij is the average j- 
hole effective mass. In this way, the contribution of each 
hole band to the polaron energy is taken into account 
explicitly in the multiband calculation. Let us notice 
that the resulting polaron correction is relatively small 
(not greater than 8%) in all cases, as shown in the tables 
below, despite the polar nature of these materials. This 
is presumably due to the fact that the coupling constant 
associated to. each hole band is relatively small (< 1.5) 
in all cases. t3 



VI. RESULTS AND DISCUSSION 

Since the reported values of effective mass parameters 
obtained by different approaches for both the RSP and 
LK Hamiltonians may have significant discrepancies, we 
have used different sets of parameterizMkms in order to 
compare the resulting impurity states.cJ~E3 

For_the wurtzite system (Tables II and III), we use Kim 
et a/.Ed RSP parameterizations obtained by full potential 
linearized muffin-tin orbital (FP-LMTO) band structure 
calculations, in which the spin-orbit coupling effects were 
obtained via the atomic-sphere approximation. We have 
also used the Suzuki et alE3 RSP parameters obtained 
by full potential linearized augmented plane wave calcu- 
lations. (FLAPW) ; a different set reported by Majewski 
et alLJ based on the norm-conserving pseudopotential 
plane-wave (P-EPW) method, and a fourth set obtained 
by Yeo et al.,c3 who employ an empirical pseudopoten- 
tial method. Notice that differences in parameters be- 
tween these two groups are typically small, but can be 
substantial in some cases (such as the value of the crys- 
tal field splitting Ai), having important consequences on 
the binding energy calculations, as we see later. 

In the case of zincblende structures, the LK— hole- 
parameters used are those reported by Kim et al.fB and 
Suzuki et a^.pl as mentioned above, and a third set by 
Wang et al.pB based on pseudopotential calculations. 
These parameters are summarized in Table IV. 

We first examine our results for the acceptor levels in 
WZ nitrides. We should emphasize here that the ex- 
perimental values of the acceptor levels in WZ-GaN are 
not without controversy. Nevertheless, in order to have 
a trend of the binding energies for different dopants, we 
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compare our theoretical calculations with the experimen- 
tal values in the literature. For GaN the results are 
listed in Table V (theoretical binding energy values are 
reported here to the nearest meV, but are calculated with 
much higher numerical accuracy for each set of param- 
eters). We note that in general the binding energies for 
different impurities are in good agreement with those val- 
ues observed in experiments-. For instance, our calcula- 
tions with the Suzuki et alrB parameters give rise to a 
binding energy for Bee an d MgG Q (241 and 253 meV, 
respectively, with the polaron correction included) which 
would seem to be in better accord with the reported est* 
perimental value (250 meV). Indeed, Salvador et alM 
reported recently room-temperature photoluminescence 
spectra of Be-doped GaN films. They found strong fea- 
tures in the 390-420 nm range which were attributed to 
the acceptor state formed by Be at about 250 meV above 
the valence band edge. Even though residual impuri- 
ties could be responsible as well for this level, no experi- 
ments have been reported_to confirm either claim. Very 
recently, Bernardini et clIlB using first principles calcula- 
tions, predict that Be is a shallow acceptor in GaN with 
a binding energy (BE) of only 60 meV, in clear contrast 
with our calculations and with the experimental data. It 
is interesting to note, however, that our BE's for Mg (~ 
200-250 meV) are in satisfactory agreement with those 
theoretical values obtained from first principles calcula- 
tionSpby Fiorentini et al&3 (~ 230 meV) and Neugebauer 
et alH (~ 200 meV). 

In contrast, the binding energies with Ref. |26f| param- 
eters, for Zn and C impurities, are overestimated with re- 
spect to the experimental values (presumably due to the 
high value of crystal field reported in p6fl). In principle 
we should expect the best fit precisely for these impuri- 
ties since they are isocoric with Ga and N, respectively, 
which would produce negligible local relaxations and core 
polarization effects. The best agreement occurs when we 
use the parameters from Kim et al.fB suggesting that 
their parameter set is somewhat better. For example, for 
Zn^a in GaN, we obtain a BE of 331 meV using Kim's 
parameters, which is in good agreement with the exper- 
imental value of 340 meV, and in excellent agreement 
with the theoretical value reported by Bernardini et al. 
(330 meV).ll3 Concerning the Cn substitutional impurity 
in a N site, we find that with exception of Suzuki's pa- 
rameters, all the hole-band sets give BE's (223-240 meV) 
comparable with the experimental value of 230 meV from 
Fischer et alLn Note that using Kim's parameters gives 
the acceptor level just even with the experimental value, 
in a nice but probably fortuitous agreement, considering 
the possible sources of systematic errors. Boguslawski 
et alu had predicted also an ionization energy for Cn 
of ~ 200 meV, while Fiorentini et alS3 report a deeper 
(~ 600 meV) value. The formation energy for this impu- 
rity is found to have also a substantial difference (1.4 eV) 
between those authors. The relatively higher relaxation 
effects predicted by Ref. || seem to play a more cru- 
cial role here. Similar discrepancies are found between 



the present work and other calculations for C&Ga and 
Sijv-u We found that C&Qa has its acceptor level (~ 260 
meV), close to that of the Mg. It is interesting to no- 
tice that temperature-dependent Hall measurements of 
Ca-doped GaN have shown that the thermal ionization 
energy level of Ca X&j_G.17 eV) is similar to that found 
in Mg (~ 0.16 eV)£Zlt_l This could indicate that the ac- 
ceptor binding energy for Ca is also close to that for Mg 
as we have indeed predicted. Similarly, Sijv was found to 
have a rather shallow level in WZ-GaN at about 0.2 eV. 
While the donor behavior of Si is well known, no reliable 
experimental evidency of Si acceptor has been reported. 

The collection of results discussedr-above indicates that 
the parameterization of Kim et al£B leads to acceptor 
binding energies in overall better agreement with the ex- 
periments and other theoretical estimates. Notice how- 
ever that the differences in binding energies in GaN with 
other sets of parameters are not large in most cases, 
within a few percent from each other. 

We would like now to comment on the effect of the 
crystal field splitting on our calculations. Whereas recent 
experiments-seems to indicate that the Ai value is about 
10 meVj23~E3 the theoretical estimates are still contro- 
versial, varying betweeri|-22-23 meV for GaN depending 
upon the approach used.c3~c3 For example, Ref. j25l and 
Pq had obtained Ai =36, and 73 meV, respectively. The 
former authors attribute the large theoretical discrepancy 
to the use of an ideal-cell internal structure parameter u 
in Ref. pq] , instead of the relaxed one. In any case, to 
illustrate the effect of the binding energies upon the Ai 
value, we show in Fig. 1 their dependence on this param- 
eter for Mga a in WZ-GaN, over a wide range. A rather 
monotonic behavior is seen in the binding energies, as 
one would expect. Note that for all parameter sets (with 
exception of those in Ref. ^7j) the BE's are consistently 
close for each Ai value. Note that using the experimental 
value of 10 meV for Ai would produce smaller binding 
energies, giving values of about 0.19 eV, regardless of the 
set employed. The behavior for other dopants shows an 
analogous trend, where the energy shift on the binding 
energy is nearly the difference in Ai values. This discus- 
sion indicates that additional experimental evidence for a 
smaller Ai value, and comparison with better optimized 
estimates, would be of interest. 

The results for A1N in the wurtzite structure are given 
in Table VI. The first thing to notice here is that due per- 
haps to the large discrepancy in Ai .scalues, —215, j=?19, 
and —58 meV for Kim,E3 Majewski,E3 and SuzukiJ^S re- 
spectively, the binding energies differ by almost a factor 
of two for different parameter sets. Notice further that 
values of ^5 and A§ also differ substantially for different 
authors, strongly affecting the band mixing and corre- 
sponding binding energies. Given the better agreement 
of Kim et al. parameters in WZ-GaN, we are inclined 
to think that the corresponding results in WZ-A1N will 
be perhaps closer to the experimental results. Unfortu- 
nately, as we mentioned earlier, the experimental spec- 
trum for acceptors in A1N is unknown at present (due 
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to the well known difficulties in doping this materialu). 
Further scrutiny of the parameters reported by these and 
future authors should be carried out to solve the disagree- 
ments. Notice that the BE of C N in WZ-A1N is found 
to exceed 0.65 eV in our calculations for all three sets of 
parameters (not shown in Table VI). This value, in the 
limit of validity perhaps of our EMT calculations, sug- 
gests nevertheless that such impurity will yield a some- 
what deeper level than those reported in Table VI. Al- 
though substitutional impurity calculations do not rep- 
resent a strict test of the band parameterizations, the 
subtle interplay of the different valence bands on the re- 
sulting binding energies (or even excited impurity states) 
provide an interesting overall consistency check. 

For the ZB phase, we notice that predicted binding 
energies are consistently smaller (by nearly a factor of 
two) than in the WZ structure of GaN. Indeed, typical 
differences of roughly 100 meV are found in the binding 
energies between the two phases (ZB and WZ) in this 
material. This would have important consequences in 
electronic uses once doping of ZB phases is stabilized. 
Concerning the resulting impurity binding energies for 
GaN, we observe that the LK parameters given by Ref. 
p5|, p{|, and (29| give rise to binding energies which 
are in close agreement with each other. We should also 
comment that a different set of band parameters in the 
ZB phase has been given by Mepey et al, using a semi- 
empirical perturbative approach.E3 However, using these 
parameters result in BE's much smaller than those pre- 
sented here. This difference, even greater in the bind- 
ing energies for ZB-A1N, reflects the more approximate 
nature of the parameters in Ref. pCf] . Notice that the 
Luttinger 7-parameters in Kim et al. are slightly smaller 
than for Suzuki et al. (or equivalently, slightly larger ef- 
fective masses) , which would be expected to yield slightly 
larger BE's for the former set of parameters, as is clearly 
seen in Table VII. ■— ■ 

Recent PL spectra of cubic GaN by As et alz3 had 
claimed as indeed we have predicted in our calculations, 
that acceptor BE's for cubic-GaN may have energies shal- 
lower than these in wurtzite-GaN. Acceptor energies of 
about 130 meV have been estimated by them. This in 
very good agreement with our calculations; as we can 
see in Table VII, the BE's ranges from ~ 130 meV for 
Si to ~ 180 meV for Zn. This acceptor level has not 
been identified and it is probably produced by residual 
impurities. 

The smaller binding energies in ZB, with respect to im- 
purities in the WZ structure, is an interesting result that 
should be understood in terms of the different band struc- 
ture parameters. Notice, however, that the difference in 
the effective Rydberg energy for WZ and ZB GaN, is not 
large at all, as seen in Tables II and IV. Similarly, the ef- 
fective Bohr radius for both structures is nearly the same, 
as illustrated in the fact that 71 is of the same order in 
both cases, and that the dielectric constant for both poly- 
types has been taken as e Q = 9.5. The polaron correction 
is certainly relatively small also, and is therefore not a 



possible source of the binding energy difference in these 
polytypes. However, the parameter that apparently gives 
rise to these large shifts in the acceptor energies could be 
identified with the in-plane heavy hole mass, which is in- 
deed larger in wurtzite than in zincblcndc for both GaN 
and A1N, and hence produces larger binding energies. In 
order to verify the effect of the different effective masses 
in the two polytypes, we have calculated the acceptor 
levels for WZ-GaN using the quasicubic sets of parame- 
ters of Ref. (25) , with the same crystal field splitting than 
the obtained for the non-quasicubic set. It turns out that 
the binding energies are smaller correspondingly, which 
confirms our assumption. One should also mention that 
just as seen in Fig. 1, a vanishingly small Ai (as is the 
case in ZB) would produce an even smaller binding en- 
ergy for a given impurity. [This would also explain the 
agreement among the three sets of parameters, since Ai 
differences are the most significant for different authors.] 
We then conclude that it is in fact a combination of the 
crystal field splitting and slightly larger hole masses that 
produce larger binding energies in WZ than in the ZB 
structure. An interesting and important effect of the dif- 
ferent lattice and band structure. 



VII. CONCLUSION 

We have carried out calculations for the shallow ac- 
ceptor energies associated to different substitutional im- 
purity atoms in GaN and A1N hosts. The calculations 
were performed within the effective mass theory, taking 
into consideration the appropriate valence band Hamil- 
tonian symmetries for the WZ and ZB polytypes, and 
using the full 6x6 acceptor Hamiltonian and included 
the actual spin-orbit energy splitting. In addition, the 
impurity pseudopotential and the electron-phonon (po- 
laron) correction has been explicitly considered. These 
more realistic treatment allows us to compare directly 
with the observed data and verify that our calculation 
produces the appropriate 'chemical shifts'. Indeed, our 
calculations of the acceptor binding energies are in quite 
good agreement with PL experiments, as the introduc- 
tion of the impurity pseudopotential seems to be an ex- 
cellent model to describe the chemical shifts associated 
with each impurity atom. It is interesting that the good 
fits were found without any adjustable parameters in the 
calculation, once the contribution due to the electron- 
phonon polar interaction was included. We find that 
small differences in the hole effective mass parameters 
could lead to relatively large discrepancies in the binding 
energies. Our overall evaluation of parameters suggests 
that the better BE values are obtained with those in Ref. 
p5[ . Correspondingly, we refer the reader to the first line 
in each impurity case in Tables V, VI, and VII, for what 
we consider the best BE estimates, within a few percent 
error. Further refinement of experimental values would 
be desirable to set narrower constraints on the theoret- 
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ical values. We also find that the binding energies for 
acceptors in the ZB structures are much shallower than 
their counterparts in the WZ structures, suggesting per- 
haps much more efficient carrier doping in those systems 
(yet to be observed experimentally). 

Finally, we should mention that preliminary studies of 
the strain effects on the acceptor binding energies show 
an increase as the strain increases, although with a much 
stronger dependence than in other III-V materials. A 
complete report of these studies will be presented else- 
where. 
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C Mg Zn 



Exp. 27.0 28.7 30.6 

Theor. 27.4 27.7 28.3 



TABLE I. Comparison between the experimental and the 
EMT-pseudopotential model of the acceptor binding energy 
for various impurities species in GaAs (e — 12.4). The en- 
ergies are in meV. The experimental values and the band 
parameters are taken from Ref. [MOj. 



WZ-GaN 



Ref. 


[§ 


n 


[§ 


{§ 


At 


-6.4 


-6.27 


-6.4 


-7.24 


A 2 


-0.5 


-0.96 


-0.8 


-0.51 


A, 


5.9 


5.70 


5.93 


6.73 


A 4 


-2.55 


-2.84 


-1.96 


-3.36 


A 5 


-2.56 


-3.18 


-2.32 


3.35 


A, 


-3.06 


-4.96 


-3.02 


-4.72 


Ai 


36 


73 


24 


22 


A 2 


5.0 


5.4 


5.4 


11/3 


A 3 


5.9 


5.4 


6.8 


11/3 


7i 


2.91 


3.80 


2.76 


3.87 


K 


51.8 


39.7 


54.6 


39.0 



TABLE II. The Rashba-Sheka-Pikus valence band param- 
eters for wurtzite GaN. The hole parameters Ai are in units 
of Tl 2 /2m , while 71 is dimensionless; Ai values represent the 
energy splittings in meV; E* is the effective Rydberg energy 
in meV. We use e = 9.5 as the dielectric constant in GaN. 
Signs of A$ and A% parameters of Ref. ]2f| have been changed 
to be consistent with those in the definition of the usual RSP 
Hamiltonian. 



WZ-A1N 


Ref. 




@ 


@ 


@ 


A ± 




-3.86 


-4.06 


-3.82 


A 2 




-0.25 


-0.26 


-0.22 


A:, 




3.58 


3.78 


3.54 


A 4 




-1.32 


-1.86 


-1.16 


Ar, 




-1.47 


-2.02 


-1.33 


A 6 




-1.64 


-3.04 


1.25 


Ai 




-215 


-58 


-219 


A 2 




6.8 


6.8 


6.6 


A 3 




5.7 


6.8 


6.7 


7i 




1.57 


2.12 


1.38 


K 




119.9 


88.8 


136.5 


TABLE III 


. The Rashba-Sheka-Pikus valence band param- 


eters for wurtzite A1N. We 


use e D = 8.5 as 


the dielectric con- 


stant in A1N. 


Parameters have the same units as indicated in 


Table II. 


Notice the enormous discrepancy in the crystal field 


splitting Ai between Ref. 


|25j, |7f, and | 


6|. 






ZB-GaN 


ZB-A1N 


Ref. 


{§ 


@ 


H 11 


{§ 


7i 


2.46 


2.70 


2.94 1.40 


1.50 


72 


0.65 


0.76 


0.89 0.35 


0.39 


73 


0.98 


1.07 


1.25 0.59 


0.62 


A D 


19 


20 


17 19 


20 


K 


51.3 


55.8 


49.3 134.5 


125.5 


TABLE IV. 


The Luttinj 


^er-Kohn valence band parameters 


for zincblende GaN and A1N. Here the dimensionless 7, are 


the hole band parameters; A Q is the enerj 


jy splitting due to 


spin-orbit interaction at the P point, and _B* are given in meV. 







WZ-GaN 


Impurity site 


e; 


E? 


E b (Exp.) 


Be Ga 


193 


204 






233 


241 


250 a 




195 


208 






185 


193 




MgGa 


204 


215 






245 


253 


250 6 




208 


221 






197 


204 




Zn Ga 


321 


331 






411 


419 


340 6 




394 


406 






352 


360 




C&Ga 


248 


259 






297 


305 






264 


276 






247 


255 




Cjv 


220 


230 






264 


272 


230 c 




228 


240 






214 


223 




Siiv 


192 


203 






231 


239 






193 


205 






183 


191 





WZ-A1N 


Impurity site 


K 


Ef 


Be^i 


223 


262 




446 


472 




283 


253 


MgAl 


465 


514 




758 


795 




721 


789 




219 


255 




438 


464 








C&Al 


204 


240 




376 


402 




203 


273 


Sijv 


214 


250 




415 


441 




245 


315 



TABLE VI. Calculated acceptor binding energies for differ- 
ent impurities in wurtzite A1N. Binding energies are ordered 
in descending order for parameters from Ref. Jjsj ], J2(|, and 
[ B7| , respectively. The large discrepancy in the calculated val- 
ues is mostly due to the important differences in the crystal 
field splitting used. 



a Reference [31] 
^Reference [4] 
Reference [32] 

TABLE V. Comparison between the calculated acceptor 
binding energies and experimental values for different substi- 
tutional impurities in wurtzite GaN. E£* and El denotes the 
estimated binding energies with and without the polaron cor- 
rection. All energies are in meV. The binding energies are 
obtained with the band-parameters from Ref. |2^], |26| , [ p7| , 
and [Q, respectively, arranged in descending order for each 
impurity. 
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ZB-GaN 



ZB-A1N 



Impurity 


El 


E* h * 


K 


E? 


Be 


124 


133 


265 


292 




117 


125 


248 


273 




126 


133 






Mg 


130 


139 


333 


360 




123 


130 


305 


330 




133 


140 






Zn 


170 


178 


261 


288 




155 


162 


245 


269 




177 


184 






Ca 


153 


162 


242 


268 




1 4^ 




227 


9^9 




157 


164 






C 


138 


147 


353 


380 




130 


138 


320 


345 




141 


148 






Si 


123 


132 


255 


281 




117 


125 


239 


264 




125 


132 







TABLE VII. Acceptor states for zincblende GaN and A1N. 
The three values shown for impurities in GaN correspond to 
those calculated with the parameterizations given by Ref . |^5| , 
[ po| , and | fp| , respectively. The two values for A1N correspond 
to Ref. Eland §1. 
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FIG. 1. Binding energies vs. crystal field splitting Ai, in 
the MgGa WZ-GaN system for different parameterizations. 
Different symbols correspond to binding energy values ob- 
tained from the effective mass parameters, o for Ref. Ef|, A 
for Ref. @, □ for Ref. 0, and O for Ref. @, respectively. 



11 



